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Abstract 

Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft 
jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin & 
A. Ajdari ( Phys. Rev. Lett. 103 , 036001 (2009)), we show that cooperativity effects resulting from the 
non-local nature of the fluidity (inverse viscosity), are intimately related to the emergence of shear-banding 
configurations. This connection materializes through the onset of inhomogeneous compact solutions (com- 
pactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. 
Compactons coexistence with regions of zero fluidity (“non-flowing vacuum”) is shown to be stabilized by 
the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity 
field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material. 
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I. INTRODUCTION 


Soft amorphous materials, including emulsions, foams, microgels, and others, display complex 
flow properties, intermediate between the solid and the liquid state of matter [1-3]. The response 
of such systems to an external shear stress shows that beyond a critical yield stress, oy, they flow 
like liquids whereas, below oy, they remain jammed and respond elastically. These materials 
are commonly denoted as yield-stress fluids (YSF). The yielding behavior makes such systems 


both interesting for applications anc 
equilibrium statistical mechanics [1 


challenging from the fundamental point of view of out-of- 
3]. Experimental studies based on the steady state rheology, 


and the behaviour in the proximity of the yield stress, have recently pointed out the existence of 
two apparently distinct classes of YSF’s. The first class |4-^], usually associated with underlying 
repulsive interactions [70, is characterized by a continuous and smooth transition between solid¬ 
like and liquid-like behaviour, as the stress is increased above ay. The corresponding rheological 
response is usually well fitted by generalised Herschel-Bulkley (HB) flow-curves, corresponding 
to a monotonous trend of the mean value of the stress a as a function of the imposed shear S : 

o = o y +AS n . (1) 

In the above equation, A is a plastic viscosity and n is a non-universal scaling exponent generally 


below 1, and often close to 0.5 [8- 


10]. Under the condition of a constant imposed shear, the steady 


flowing state of these materials, often denoted as “simple” YSF, is homogeneous in space, even 
at vanishingly small imposed shear. The second class of YSF’s displays a discontinuous transi¬ 


tion from so 
threshold [7, 


id-like to liquid-like behaviour upon increasing an imposed stress above the critical 
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14]: the shear jumps discontinuously from a zero value below the threshold stress 


to a finite value just above it. This may also be regarded as a discontinuous jump in the viscosity, 
from formally infinity, just below threshold, to a finite value just above it, a phenomenon often in- 
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150 . This scenario underlies the formation 


dicated by using the term viscosity bifurcation [3J, 4, 
of shear-bands, where fluidized zones can coexist with solid-like regions, under the conditions of a 
constant imposed shear. These kind of materials are usually referred to as ’’viscosity-bifurcating” 
YSF’s. Within the above two classes of materials, only viscosity-bifurcating YSF’s appear to sup¬ 
port shear-banding as a permanent response to a steady imposed shear. 

Developing predictive theories for the deformation and flow of amorphous materials, study the 
formation of shear-bands, as well as identifying the correct structural variables for these systems 
remains an ongoing challenge: in the absence of a comprehensive microscopic theory, various 
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mesoscopic models have been vigorously pursued in the literature [ 16f-28]. It is also commonly 
accepted that as soon as the flow becomes heterogeneous, a description of the rheological be¬ 


haviour solely in terms of the flow-curve is insufficient [28 


-36]. This has been illustrated in the 


experimental work of Goyon et al. [29D, who showed that a single flow-curve is not able to account 
for the flow profile of a concentrated emulsion in a microfluidic channel. That finding triggered 
the need to properly bridge between local and global rheology of the soft-glassy. Specifically, 
Goyon et al. [29[] introduced the concept of spatial cooperativity length !q c , by postulating that 
the fluidity, f — S/a, defined as the ratio between shear S and shear stress a. follows a non-local 
diffusion-relaxation equation when it deviates from its bulk value 


+fb(o(r)) -fir) = 0. 


( 2 ) 


The quantity fj, = fif e (7) ) is the bulk fluidity, i.e., the value of the fluidity in the absence of spatial 
heterogeneities. The bulk fluidity depends upon the local stress < 7 , whereas / = fir) depends upon 
the position in space. Its value is equal to fb without the effect of cooperativity (q c — 0). The 
spatial cooperativity tq c has been shown to be in the order of few times the size of the elementary 
microstructural constituents flEjlQQH. The non-local equation © has been justified fl] 
by using a kinetic model for the elastoplastic (KEP) dynamics of a jammed material, which takes 
the form of a nonlocal kinetic equation for the stress distribution function. Such model predicts 
nonlocal equations of the form ©, plus an equation predicting a proportionality between the 
fluidity and the rate of plastic events F, 


/ = - oc r. 

a 


(3) 


An interesting interpretation of the diffusion equation © has been put forward in [30], based on a 
“free-energy” for the rate of plastic events T 


f{ r)= / ( ©(r^j + ^-ivri 2 )d 3 x 


(4) 

where the bulk potential <u(r, a) embeds the information about the non-linear rheology, whereas 
an inhomogeneity parameter q appears as a multiplicative factor in front of the gradient terms 
and may be directly related to the cooperativity length tq c . Upon minimization, the free-energy 
formulation © leads to the HB form of a flow-curve © in case of homogeneous flow, and to the 
fluidity equation © in case of non-homogeneous solution. As pointed out in [30], the square- 
gradient expression of the free-energy © opens up a far-reaching connection to shear-banding. 
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In this picture, shear-banding would correspond to a first order phase transition scenario: i.e., the 
spatial coexistence of two states of different fluidity values for the same shear stress. It is the 
purpose of the present paper to elaborate on these concepts, proposing a new viewpoint on the 
formation of shear-bands in soft matter systems. Shear-banding will be linked to the onset of 
compact configurations of the fluidity field (compactons), which correspond to the local minima 
of the free-energy functional. More specifically, we start from the non-local formulation due to 
Bocquet et al. [30], whose main idea is to introduce the fluidity as the order parameter of a corre¬ 
sponding free-energy. Free-energy minimization leads to a non-linear Helmoltz equation, whose 
solutions describe spatial relaxation to a uniform background fluidity, corresponding to homoge¬ 
neous bulk rheology, typically in HB form. We also inspect the dynamics of the system subject 
to stochastic perturbations, arguably related to a form of mechanical noise in the system. None 
of these concepts is brand new in the literature, although the deep and non trivial consequences 
shown in this paper definitely are. We are going to show that the cooperative effects, resulting 
from the non-local nature of the fluidity, are intimately related to the emergence of shear-banding 
configurations. This connection materializes through the onset of inhomogeneous compact solu¬ 
tions, wherein the fluidity is confined to finite-support sub-regions of the flow and zero elsewhere 
(non-flowing vacuum). In the absence of noise, compactons attain lower free-energy minima than 
the homogeneous HB solution, thereby realizing metastable shear-bands. Indeed, since the non¬ 
flowing region is unstable, such shear-bands cannot survive indefinitely: depending on the initial 
conditions, shear-bands may or may not occur, but even when they do, in the time-asymptotic limit 
they surrender to homogeneous HB configurations. The picture takes a drastic upturn once noise 
is taken into account. Here, a new qualitative effect arises: the effective free-energy, including 
renormalized fluctuations, develops two local minima, corresponding to a stable coexistence of 
compactons and non-flowing vacua. Under such conditions, permanent shear-banding solutions 
can indeed be observed. The emerging picture is conceptually sound and appealing: compactons 
represent natural carriers of shear-bands. They attain local minimization of the free-energy func¬ 
tional, but in the absence of noise they ultimately surrender to homogeneous configurations due to 
the instability of non-flowing vacuum. 

We remark that the ’’minimization of the free energy” is a (quite) strong assumption that can be 
hardly justified at this stage. Nevertheless, we shall show that this idea has far reaching conse¬ 
quences. In particular, as a consequence of the noise, there is a renormalization of the free energy 
leading to an unstable rheological branch. Interestingly, the renormalization effects can be com- 
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puted analytically within the framework of a Hartree-like approximation (see Appendix®. The 
analysis can be performed for a very simple structure of the noise, namely a white noise in space 
and time. We may argue that such a choice of the noise can be considered quite unrealistic and 
in principle the noise could depend on a external forcing or the rate of energy dissipation and it 
may be correlated in space and/or in time. Our choice has been done to maintain the whole theory 
analytically manageable without introducing a larger set of independent parameters. Based on 
our results, future investigations may include more realistic structures of the noise term and the 
relative renormalization effects. 

The the paper is organized as follows: in Section [TT] we describe the basic features of the free- 
energy model, specializing to the geometry of the Couette flow; in Section [ITT] we explore the 
energy landscape of the model and we describe the compact solutions; in Section [TV] we study the 
free-energy of both the HB and compact solutions and the resulting velocity profiles; the stabil¬ 
ity of HB and compact solutions, as well as the aging properties of the model, are the subject of 
Section |V] In Section |VI] we analyze the behaviour of the system under stochastic perturbations 
(noise). The far-reaching implications of the stochastic perturbations for the geometry of the Cou¬ 
ette flow are illustrated in Section lYlIl In Section IVIIIl we offer a discussion with comparisons 
between our model and other earlier works that report about the formation of permanent shear- 
bands. Conclusions follow in Section [[X] Technical details for the stability of compact solutions 
and the self-consistent field approximation in presence of noise are reported in Appendices lAllBl 


II. MODEL EQUATIONS 


We are interested in the dynamics of a soft-glassy system within two walls at y — 0 and y — L 
driven with a velocity difference A U = SL at the boundaries, with S the imposed shear. Using a HB 
relation with n — 1/2, we rescale the original variables according to the following transformations: 


y~> 


y 




o ->■ 


L \ A / oy 

As a consequence, the HB relation ([[]) retains a unitary yield stress 


(5) 


(7=1 +S 1//2 . 


( 6 ) 
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Following Bocquet et al. [30|], we write equation © as 


F \f\=[ (“(r,®) + - 


ivn' 


dy — 


r=/ 


where 


(C7-1) 


(7) 


dy 


m(o) = 


rV2 


( 8 ) 


Equations ([7])-® are the starting point of our investigations. In case of a spatially homogeneous 
solution, the minimum of F[f] is given by 

(o- l ) 2 


f = m = 


(9) 


that is the FIB solution. For spatially non-homogeneous solutions, upon linearizing around fo = 
nr, we get the equation 

gdyyf-f + f b = 0 (10) 

which directly maps into ©with a squared cooperativity length t, 2 = 2c, 2 /m diverging at the yield 
stress (as a —> 1, m 0) [30]. Let us remark that the value of m — (<J — 1)/(7 172 in the above 
equations is consistent with the definition of the fluidity f — S/o. In a Couette flow, assuming as 
usual that the stress is constant in space, we can obtain a constraint between the imposed shear and 
the space-averaged fluidity 


S = a f(y) dy. 


( 11 ) 


III. ENERGY LANDSCAPE AT IMPOSED STRESS: HERSCHEL-BULKLEY VS. COMPACT 
SOLUTIONS 


Let us start with the analysis of the energy landscape in the system at an imposed constant stress 
a. We want to specialize our analysis to those situations where the “order parameter” / is positive 
definite. This request can be ensured if we set f — (f> 2 and use the corresponding free-energy 
functional, directly obtained from © with the substitution / = (j) 2 : 

r 1 


F[$\= \ L(4>,d y 4>)dy = 




dy. 


( 12 ) 
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Note that the term / 5 / 2 has been rewritten as (j) 4 \(j)\ to guarantee F[<j>] > 0 in the limit (j) —y ±°°. 
We then look for the local extrema of the free-energy functional (fl2l) . The variational equation 

If = 0 § ives: 

2E, 2 <\> 2 dyyty +2^ 2 <^(5y<^) 2 + m(cr)0 3 — 0 3 |0| = 0. (13) 

Equation (fl3l) exhibits solutions with constant order parameter (</),„ = m ), i.e., the bulk HB solu¬ 
tions. However, also compact solutions are possible 


0c (y) 


0o(y) y < > 3^i 

0£<» yo < >’ < >'l 


(14) 


with non zero values of the order parameter only in some compact sub-domain, say [yo,yi], an d 
zero elsewhere (see the top panel of Figure Q] for a sketch). In the above, 0o(y) = 0 is a “vacuum” 
field, corresponding to a zero-fluidity (non-flowing) state. The structure of the compact solution 
in the region where the order parameter (j) is different from zero, i.e., the function 0^(y) (hereafter 
named compacton ), can be further characterized by taking one quadrature of equation (fl3l) 

<§ 2 0 2 (<5>'0) 2 = E- ^^-0 4 + ^0 5 (15) 


where E is a positive constant and where, for the sake of simplicity, we assume (j) > 0. By increas¬ 
ing 0 from zero to larger values, and for sufficiently small values of E [ 84j] 

1 (<7-l) 5 


E <E m (o) = 


(16) 


20 a 5 / 2 

the r.h.s. of (fl5l) first becomes zero for a value of the order parameter denoted with (j) — 0 < m. 
Then, we can look for a solution of (fl5l) localized in the interval [yo,yi], with 0(y) = 0 and 
y — (yo + yi)/2. Based on the structure of the exact solution 


V 


io 


E-'EiEly 4 + 1^5 


d\ff = 


y-y o 


(17) 


4 ' 5 

we obtain the characteristic size of the region, |yi — yo I = 2 l c , where the order parameter is different 
from zero 


k = 4 / dy 


¥ 


E-'EiEly 4+Il/r 5 


(18) 


4 V 1 5 

It is easy to check that for y close to yo, i.e., for small (j), we have 


0 


1/4 

^ ^ (y-yo) 1/2 


r\j 


E 

V 
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(19) 












implying that d y (j) diverges as I /(y — yo) 1/ 2 at small (j). Importantly, equation (fl3l) is well defined 
even at the singular points since all the divergences cancel out. It is also important to highlight 
that the characteristic size of the compact region (fl8l) scales proportionally to c ,, hence no compact 
solution can be achieved without cooperativity. Since compact solutions do not overlap, superpo¬ 
sitions of compactons still correspond to local extrema of the free-energy. This implies that the 
energy in equation (fl5l) takes the form of a piece-wise constant function, which attains distinct 
non-zero values in different compactons and is zero elsewhere. 


IV. FLOW PROFILES AT IMPOSED SHEAR: LINKING COMPACTONS TO SHEAR-BANDS 


As discussed in the previous Section, the energy landscape of our system is characterized 
by various stationary solutions corresponding to both homogeneous (HB solution) and non- 
homogeneous (compact) solutions. Now, we study the corresponding free-energy of those so¬ 
lutions and the resulting velocity profiles. To that purpose we go back to the geometry of the 
Couette flow at imposed shear and consider for simplicity the case of a compact solution with a 
compacton adjacent to the upper wall of the channel with the property </>e( y) > 0 for l c <y< l 
and d y (j) — 0 at y = 1 (see bottom panel of Figure [2]). This choice of boundary conditions makes 
it easy to discuss the case of compact solutions close to the boundaries, although it differs from 


the usua 


fluidity [2' 


choice of Dirichelet boundary conditions, where a new parameter representing the wall 


5 , 


28 


38] is introduced. In the Couette flow with imposed shear S, the overall free- 


energy of both the compact and HB solutions depends on the shear S as follows: given the stress 
a, the full non linear solution of equation (fl3T) must be found for a value of m(o) = (a — 1)/<T 1//2 
consistent with 


t7 / </> 2 (y) dy = S. 
Jo 


( 20 ) 


For the HB solution, corresponding to (j) = (f) m — m, the constraint (l20l) imposes a — S/in 2 so that 
the HB relation ([8]) implies m(o(S))= ( 1+ si/ 2 p /2 an d the corresponding free-energy is given by 

1 S 5 /2 


Fhb(S) = 


( 21 ) 


10(l+S 1 /2)5/2- 

The free-energy for the compact solution is different. Given the shear S, the free-energy is not 
uniquely determined, as multiple choices of a and E are compatible with the constraint (l20l) . This 
makes the free-energy dependent on both o and S: once these two parameters are fixed, we are 
able to determine the appropriate E in dT31) that makes possible to derive the ()/.■ that satisfies the 
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constraint (l20l) . For a given stress a, this obviously leads to a larger l c at increasing S, i.e., the 
size of the compacton is increased to verify the constraint (l20l) (see the bottom panel of Figure 
[[])• In the top panel of Figure [2] we show the free-energy Fhb(S) of the FIB solution and the free- 
energy of the compact solution F C (S, <j) as a function of the shear S. For the compact solution, we 
choose two different values of the stress, a = 1.4 and a — 1.5. We see that the cases with compact 
solutions show a free-energy smaller than the HB value up to a critical shear S cr which is stress 
dependent 

F c (S,o)<F hb (S) S<S cr (o). (22) 


In correspondence of S cr , the size of the compact region becomes of the order of the channel size. 
Compact solutions in y imply that the shear, according to the definition of the fluidity f — (j) 2 . 
is different from zero only in a compact region. This suggests, by itself, an intriguing link to 
shear-banding , since it permits coexistence of (compact) flowing and non-flowing states within 
the same spatial flow configuration. This can be evinced from the bottom panel of Figure |2l where 
we analyze a situation with imposed shear S — 0.1 and stress o = 1.5. 

Central to this free-energetic picture is the imposition of the global constraint (l20l) . which is key to 
tip the free-energy balance in favor of compact solutions versus HB solutions, F C ((J,S) < Fhb(S), 
in the proper range of the shear S, as reported in Figured We wish to observe that, had we imposed 
the stress o, we would have obtained 


F c (a) = 2 


ry\ 

ho 


-^-m(d)0 4 + ^ 4 |0|+£ 2 0 2 (<9y0) 2 dy > F hb {o) 


(23) 


meaning that the compact solution </> c has larger free-energy than the HB solution. In other words, 
by imposing the stress instead of the shear, the HB solution would always show up. 


V. STABILITY OF COMPACT SOLUTIONS: TIME DYNAMICS AND AGING 


The connection between shear-bandings and the compact solutions (j) c is extremely appealing, 
but we hasten to remark important features about the stability of the shear-banding solutions thus 
obtained. To this aim, we resort to the very simple equation 


d t (l) = 


8F_ 

8(f) 


(24) 


to define a time-evolution of the order parameter that is consistent with the attainment of local 
extrema of the free-energy functional in the long-time. More sophisticated formulations of the 
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dynamics may eventually be considered. 

As a matter of fact, much of the insight conveyed by Figure [2] can also be gained by analyzing the 
scaling properties of the model. In particular, equations of motion (1241) are invariant upon the scale 
transformation: 

y^A fl y £->• X b l 2 + a E, t ^ X~ 3b t m^X b m. (25) 

By using (1251) we determine the scaling properties of the free-energy as 

F ->• X 5b+a F. (26) 

The choice a — 0 and b — 1, corresponding to a linear velocity profile with increasing amplitude 
m and no shape change, implies a free-energy that scales as F ~ A//#A 5 . The choice a — 1 and 
b = 0 corresponds to compact solutions, namely increasing size of the shear-band with no increase 
in amplitude m, and delivers F ~ A C X. Note that both constants A HB and A c are negative, so that 
the minimum corresponds to the larger absolute value. Finally, under the assumption that S is 
sufficiently small, i.e., <7 ~ 1, equation (l20l) delivers Shb A 2 for a — 0 and b — 1 and S c ~ A for 
a—l and b — 0. Based on the above scaling relations, we obtain: 

Fhb~A hb S 5 / 2 (27) 

F c ~ A C S . (28) 

The result is that, up to a critical value S = S cr = 0 A c /Ahb ) 2/3 , one has F c < Fhb, so that the 
compactons are favored with respect to the homogeneous configurations. Scaling laws (l27l)-(l28T) 
correspond to the functions reported in Figure O 

Given the time dynamics (1241) . one can verify that the compactons (f>E are stable against pertur¬ 
bations of the order parameter. In particular, one can compute 8F = F[(j>E + 8^] — F [</)£■] up to 
the second order term in 8(j), and show that 8F is positive defined for an energy E sufficiently 
close to /n 5 (c7)/20 (see Appendix lAl). However, the non-flowing state 0o in (fl4l) is unstable, and 
the overall stability of the system must take into account both compact and vacuum components. 
The compacton is thus expected to grow and “eat up” the unstable non-flowing state By this 
process, a critical situation is attained whenever the compacton hits the size of the channel: at that 
point, the HB solution is energetically favored again and the compacton yields to an extended HB 
profile, going back to a situation with a global bulk rheology. Under such conditions, and consis¬ 
tently with the fact that compact regions evolve at a faster rate than the “vacuum”, it is apparent 
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that the time-relaxation of the overall compact+vacuum system becomes heterogeneous not only 
in space but also in time, i.e., the systems shows aging. To highlight this effect we have conducted 
numerical simulations: the interval [0,1] has been discretized with 512 collocation points and an 
Euler-Cauchy scheme with integration step dt = 10 -3 has been used for the time dynamics of the 
free-energy. The parameter q has been fixed to q — 0.04. At time t — 0, we start with an initial 
condition 0 = 0.2 in the region [0 : 0.8] and 0 = 0.5 in [0.8 : 1], i.e., we start with a non-uniform 
initial fluidity. We integrate the equation of motion for a time f w and, at t = / w we apply a shear 
S = 0.04 and we compute the resulting stress <r(f,f w ). In Figure [3] we show c(t,r w ): after an 
almost linear growth, the stress reaches a maximum Om(^w) and then eventually decays to the HB 


value 1 -t-S 1 / 2 at time lz,(f w )- Note that Om and tt depend on t 


le waiting time t w , in particular there 


is an overshoot that depends on the age t w of the sample [39]. Similar results have been reported 


experimentally in [5|] and in MD simulations of Lennard-Jones glasses 1140 . 


4111 (see also [42] for 


a recent review). In the inset of Figure [3[ we show ti as a function of t w . It is apparent how it is 
very well fitted by a logarithmic function of t w . We remark that the same effect does not appear if 
we consider the dynamic equation for the fluidity as 


3r/ = 


8F_ 

W 


(29) 


In fact, upon multiplying equation (1241) by 0, we can rewrite equation (l24l) in the form 

8F 


dtf = ~f 


8f 


(30) 


The comparison between equations (l29l) and (l30l) highlights the physical meaning of the relation 
f — (f) 2 : upon the assumption of a simple first order steepest-descent dynamics (l24l) . the trans¬ 
formation from / to (j) implies that the former evolves on a typical configuration-dependent time 
scale 1// (l30l) . so that non-flowing states / ~ 0 take a virtually infinite time to relax, as opposed to 
the flowing ones. This configuration-dependent time-scale separation lies at the heart of the aging 
phenomena described in Figure [3] 


VI. EFFECTS OF MECHANICAL NOISE AT IMPOSED STRESS 

In the previous Section we have shown that in a Couette flow the coupling between the external 
shear and the stress is such that compactons (j)£ are stable states and the associated free-energy 
F c (S,o ) is smaller than the HB value up to a critical shear S cr . However, the region where the 
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order parameter is zero (the non-fluidized band) is unstable, so that the formation of permanent 
shear-bands is not possible and they can be observed only for a finite (possibly long) time. So 
much for the deterministic picture. The next natural question, is to inspect the behavior of the 
system under stochastic perturbationsjnoise). We shall refrain from identifying such noise with 
any thermodynamic temperature [j 18l . I43l445i ] . Actually, we rather think of it as a mechanical noise 
due to dynamic heterogeneities. A naive expectation about the effect of the noise is that it raises 
the instability of the background field (j )o so that the system reaches the HB solution in a shorter 
time. However, this expectation is not true and we will discover something unexpected, i.e., that 
the vacuum solution - and hence the compact solutions - are stabilized by noise. For the sake of 
simplicity, we include the constraint (l20l) in the next Section whereas, here, we analyze the role of 
fluctuations at constant stress. We consider the time evolution (l24l) : 


+ a/£ w (y,t) 

F [<t>] =2/o [-\m(a)(j) 4 + ^(j) 4 \(l)\+^ 2 (f) 2 (dy(l)) 2 ]dy 


(31) 


at imposed external stress <J, i.e., at imposed m(a) = (a — 1 )/c 1 / 2 . In the above w (y,t) is a 
c)-correlated white noise in space and time 


(w(y 1 ,t 1 )w(y 2 T2)) = < 5 (yi -y2)S(h-t 2 ). 


(32) 


It must be understood that the solution of equation OTI) is defined with an ultraviolet cutoff A = 
1 /Am needed for the regularization at small scales. A first, non trivial, consequence of the time 
dynamics (OTT) is that for m — 0 (i.e., <7 < cy), the space averaged fluidity /o of the system does 
not vanish. As discussed in the Appendix |B] we have 

V £k M 


fo = <r> 


%k M 


(33) 


These results can be tested quite accurately by using numerical simulations. In Figured we show 
the value of fo as a function of E, (main figure) for £ = 10 8 , whereas in the inset of the same 
Figure we show fo as a function of £. In both cases, the scaling predicted by (1331) is extremely 
well verified supporting the theoretical results obtained in appendix [0 With the interpretation of 
the fluidity as the rate of plastic events in the system, equation (1331) tells us that a non zero value 


of the fluidity is still present at zero external forcing (i.e., S — 0 and/or t7 < cy) [4g-48]. Further 


discussions on the nature and properties of the noise will be reported in the end of the present 
Section and in the conclusions 
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Next, we discuss the effect of the noise for m > 0. Some interesting insights can be gained by 
using a self-consistent Hartree-like approximation [49H51I1. as discussed in detail in Appendix |B] 
This amounts to consider the free-energy 


F W = 2 


-i/n(<T)(f . 4 + j<t> 4 \<p\ + j<t> 2 + j(d,4) 2 


dy 


(34) 


where 


R = 2^(MY) D = 2% 2 W) 


— ie2/ A 2\ 


(35) 


need to be determined in a self-consistent way. Based on equation (HIT) and the results discussed 
in Appendix |B] we can formally write the model equations as 


( d t ty — — 24^ + 2Dd yy (j) + \/ew (>’,?) 

{v„(*) = -X<r)* 4 +}* 4 W+f* 2 . 

Owing to renormalization effects, fluctuations turn the bare free-energy into an effective one, 
whose properties may lead to qualitatively new phenomena not contained in the original formula¬ 
tion. For the case in point, the potential V flows into an effective one V eff , supporting qualitatively 
new extrema through a renormalized “mass” term f <)) 2 . Besides, a new diffusion term arises with 
no counterpart in the noise-free formulation. A full non-linear treatment would require that R and 
D were treated self-consistently, i.e., taking into account their functional dependence on the con¬ 
figurational statistics of the system. However, as we shall see, significant insights can be gained 
by provisionally treating both quantities as constant parameters, and deferring a self-consistency 
check to a subsequent numerical solution. Thus, upon assuming R and D in (l36l) as “constant” 
parameters, we look for local minima of the effective potential by solving (0 > 0): 


dVjf 

d(j> 


= (R(j> —m(j) 3 + 0 4 ) 


d<t > 2 


> 0 . 



(37) 


It is easy to show that (1371) has one solution </>„„„ — (j) 0 — 0 for m < m* ~ (277?/4) 1 / 3 , whereas 
for m > m*, two local minima appear at — (j) 0 — 0 and </>„„„ = (f> m ~ m — R/m 2 + 0((R/m 2 ) 2 ), 
together with a local maximum in (j) — (j>M- This is shown in the top panel in Figure [5] The local 
minimum 0,„ is reminiscent of the HB solution already discussed in the previous Sections: this 
state receives a normalization and its stability changes, since now it appears as a local minimum 
only for m > m*. Note that there is a range of in (m* < m < m c ) where the local minimum </>,„ 
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appears, but </>o is still the global minimum. Above m c the minimum (j) m becomes the global 
minimum and the (stable) fluidized band is energetically favored. It is also interesting to extract 
the (stationary) homogeneous flow-curves from the local extrema solutions. Beyond the vacuum 
branch (j) = <f>o, we obtain a fluidized branch from the equation R — m(f) 2 + (j) 3 = 0. In the bottom 
panel of figure [5] we show the stationary solutions, in terms of stress and shear, corresponding to 
this fluidized branch: unlike the simple HB model, in presence of fluctuations (R > 0), this curve is 
a decreasing one before it slopes upward at higher shear. By fixing m (i.e. the stress c), the system 
is able to sustain the homogeneous fluid state only for m > m*. For m < m*, the only possible 
homogeneous solution is (j) — 0q. The decreasing part of the flow-curve corresponds to the local 
maximum (j>M of V eff and therefore, in terms of our dynamical equations, it is unstable. We remark 
that we impose the stress as a spatially constant variable and this instability is different from the 
“mechanical” instability induced by inertial terms [24]. 

However, this is what we formally expect, based on the form of V eff , by treating R and D as 
constant parameters. As mentioned earlier on, the complexity arises from the fact that both R and 
D depend on the fluctuations themselves (l37l) . More importantly, R and D are state dependent , as 
can be seen by perturbative calculations (i.e., assuming £ to be small) and by properly linearizing 
the potential near the local minima. In particular, an explicit computation (see Appendix |B]) shows 
that the scaling with respect to the noise strength £ is different, depending on the local minima (0o 
or (j> m ) around which they are computed: 


~ ~ eW(2m ). 


(38) 


Consequently, the bifurcation point m = m c is different, depending on whether one linearizes 
around (j)Q or </),„. We label these two bifurcation points in i{)] and m^ l \ so defining also the corre¬ 
sponding critical stresses, <7 ( °) and cO), based on ([8]) 

<j(°) — 1 


rn^\e) — 


m 


<»(£)= 


(a(°))V 2 

— 1 


£ 1/2 £V 6 
K M 




(39) 


(40) 


(O' (1)) 1/2 M 

Due to the different scaling properties in £, we find that in 1 1 > < or equivalently that (j(i) < 
a (0) . The above picture implies hysteresis in the system in the region [cr^ 1 ), o‘(°)]. Indeed, any 
initial state close to (j)o is expected to attain the only stable solution of the system up to m = m°>, 
i.e., up to a maximum stress <7 (0) . Similarly, starting with an initial state close to (j) m , the system 
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remains close to this state only for m > nP 1 ', i.e ., only above a given stress <7 ( 1 \ Therefore, once 
the stress falls within the range a E [cr*■ a* 0 *], two stable solutions are expected, depending on the 

initial conditions. As a matter of fact, these expectations are confirmed by the results of numerical 
simulations shown in Figure [bj The system is initialized in the state (j) o and we then consider the 
time dynamics based on equation (OTT) with the noise in the range [10~ 8 ,5 x 10~ 7 ]. We slowly 
change - step by step - the stress a by increasing m and wait for the system to reach a stationary 
(in statistical sense) state, where we compute the apparent shear based on S = o(® 2 ). For this 
particular simulation, we must use as boundary condition (j) — (p w — m. The reason for this choice 
is to avoid trapping in the stable state (p — (po for any in. In the top panel of Figure [6l we report the 
stress/shear relation obtained with £ = 10 s : the red circles correspond to the simulations starting 
with the state (f >o and slowly increasing the stress; the blue squares correspond to the case where 
we start from state (p m and a relatively large value of the stress, and then decrease the stress. In the 
bottom panel of Figured we show similar results obtained with the noise amplitude £ = 5x 10 -7 . 
Both figures support a well defined hysteresis effect, which is amplified by increasing the noise 
amplitude. Note that in the region c < c ( 0, only one stable solution is attained. 

From the shear-stress curves we can also determine the corresponding values of the stresses at 
which the transition occurs. These are reported in Figure [71 where we plot both my 2 ’ (red circles) 
and mS 2 ’ (blue squares) as a function of £. The scaling predictions from equations (l39l)- (l40l) are 
also reported: note that although the two scaling laws are close to each other, it is clear that the 
numerical results cannot be fitted with the same scaling exponents in £. Our simulations therefore 
indicate that our analytical results (l39l)- (l40l) are in very good agreement with the complex non 
linear dynamics of the system. The existence of a hysteresis effect in the system is an important 
point since it is one of the many puzzling results obtained in laboratory experiments in soft glasses 


(see [7|, 


11 


14] and references therein). Here, we can state that the hysteresis effect is due to the 


noise in the system and to the renormalization effects in the dynamics. 

We can also use our findings to study other peculiar facets of the phenomenon of viscosity 
bifurcation. To this aim, we consider an initial pre-sheared state with fixed shear S a = <p 2 <J — 0.01 
at time t — 0. This is realized with a homogeneous initial condition (p (y. 0) = (p a = \/S a / (7. We 
then vary the applied stress a and study the behavior of the apparent shear S(t ) = cy / 0 ' <p 2 (y,t)dy 
as a function of time. In the top panel of Figure[8l we show S(t) for different values of the applied 
stress <7 with £ = 10 -8 , E, — 0.04. There is a clear bifurcation between the values of the stress 
below a critical value (J (c K at which the shear tends to vanish, and those above (J <C K where the 
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shear goes asymptotically to a non zero value. Similarly, in the bottom panel of Figured we show 
S(t) for different values of the initial shear S a with £ = 10~ 8 , q = 0.04 and fixed applied stress 
o — 1.2. For small values of S a , the apparent shear tends to vanish whereas for S a greater than 
some critical value, the shear goes asymptotically to a non zero value. 

Before closing this Section, we would like to remark that there is a simple argument to ob¬ 
tain the scaling laws reported in equations (l39l)-(l40l). To this aim, let us consider the probability 
distribution P[<j>] given by: 


if 1 

£ Jo 


j0 4 +yl0l ) d y 


P[(t>] = J^'exp 

It is easy to verify that (l4TT) is invariant under the scale transformation 


(41) 


$^X b/2+a $ m^rX b m £ l 5b+a £. (42) 


For a — 0 the scale transformation (1421) ignores spatial structures in the system and we obtain a 
scaling given by m ~ fi 1 / 5 . For b — 1 the scale transformation (1421) takes into account the spatial 
structure and we obtain the scaling m ~ fi 1 / 6 . Therefore the scale transformation (1421) makes it 
possible to understand why there exists two different scaling regimes (l39l)- (l40l) . It is important 
to remark that the identification of P as the probability distribution implies that the system is in 
equilibrium. 

It is worthwhile to discuss in more detail the peculiar effects of the noise in our theory. As we 
already pointed out in the Introduction, we chose a very simple form of the noise to compute ana¬ 
lytically the renormalization effects of the free energy. We remark that renormalization effects can 
be computed independently of a “free energy formulation” for the fluidity. In principle, the noise 
can depend on the external forcing and/or the rate of energy dissipation and it can be correlated in 
space and time. The effect of space/time correlations can be computed by a suitable (non trivial) 
generalization of the methods illustrated in Appendix |B] Similarly, given a prescribed formulation 
of the noise variance as a function of the forcing, it is possible to generalize our results. At the 
present stage, we wish to keep our discussion as simple as possible. We argue that the qualitative 
picture emerging from our results is quite robust and independent (qualitatively) of the structure 
of the noise. 
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VII. EFFECTS OF MECHANICAL NOISE AT IMPOSED SHEAR 


Given the findings reported in Section ED we now consider the dynamics given by equation 
(l3TI) . together with the constraint (l20l) in a Couette flow. If the vacuum 0o is stabilized by fluctua¬ 
tions, the overall stability of compact solutions (hence the co-existence of flowing and non-flowing 
stable states) discloses the possibility to support stable shear-band configurations. Based on the 
results reported in Section [VH we expect the following scenario: the vacuum branch and the stable 
fluidized branch are separated by a range of shear. By fixing the imposed shear in such a range, a 
(stationary) homogeneous solution would be unstable, as it would fall in the decreasing part of the 
flow-curve (see figure [5]). We also know that compactons </>£ correspond to a minimum of the free- 
energy, a property which stays unchanged even in presence of external noise. If so, the compact 
solution (vacuum + compacton) could then be chosen. Even in case the homogeneous fluidized 
solution becomes stable, the compacton solution would still be selected as long as it manages to 
attain a lower value of the free energy. These expectations are indeed borne out by numerical sim¬ 
ulations that highlight the way in which the stress and the size of the compacton adjust matching 
the imposed shear. Let us first consider the dynamics of the system at an imposed shear S = 0.04 
with “weak” and “sufficiently strong” noise, £ = 10“ 8 and 10~ 7 . By “strong” we imply fluidity 
fluctuations in the order of a few percent of the deterministic value / = m, whereas fluctuations 
induced by a “weak” noise are below such threshold. For each noise strength, we consider two 
initial band sizes below and above the half size of the channel, / c (0) = 0.2 and / c (0) = 0.75. Using 
(l39l) and the analysis presented in the previous Section, we estimate that a*°)(£ = 10 -8 ) ~ 1.2 
and cj ( 0 )(£ = 10~ 7 ) ~ 1.4. The numerical results obtained using the four possible combinations 
(initial conditions and noise strengths) are shown in the top panel of Figure [9l where we plot the 
time evolution of the shear-band l c {t). Let us begin by discussing the results with “strong” noise, 
£ = 10 -7 . With / c (0) = 0.75, the initial stress is below <7 ( °)(£ = 10“ 7 ) ~ 1.4. In order to mini¬ 
mize the global free-energy (see also Figure [2]), the system increases the stress, i.e., it decreases 
the size of the shear-band l c , based on (l20l) . This shrinking process goes on until l c becomes short 
enough to take o close to o l{)] . At this stage, any further narrowing of the compacton would 
take (7 beyond the threshold a (0 ), thus destabilizing the region (j> — 0. Under such conditions, the 
shear-band solution increases its size again to l c . In other words, l c is pinned down as the size 
compatible with the stress (7 = (7 (0 \ based on (l20l) . The same reasoning carries on to the case 
l c ( 0) — 0.20, in which case (7 > <7 (0 ), so that l c increases until <7 reaches the value <7 (0 ^. Either 
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ways, the only stable solution of the system corresponds to a shear-band with <y = Let us 
now turn to the “weak-noise” scenario, £ = 10 8 . In this case, no shear-band solution is expected 
because <7*°) ~ 1.2: the shear S = 0.04 would select a point on the curve F c (o — a^°\S) whose 
free-energy becomes comparable with the one of the HB solution, F c (o = <7 ( °), 0.04) « Fhb{ 0.04), 
and the corresponding homogeneous solution is stable. Indeed, from Figure (9[ we observe that, 
starting with l c (0) — 0.75, the system evolves towards metastable shear-bands, which disappear as 
soon as the compacton hits the size of the domain, after about 3000 — 4000 time units. The case 
l c (0) = 0.20 presents however a different scenario. In this case, the numerical simulations show 
that l c remains basically constant, at least in the time span covered by the simulations. The reasons 
is that, even though that band is unstable, one must wait for a longer time before reaching the state 
l c — 1. Upon decreasing Z c (0), the value of the initial stress a due to (120b increases, and so it 
does the energy barrier protecting the asymptotic state l c — 1, leading to an even longer waiting 
time before the band is destabilized. This is appreciated from the bottom panel of Figure [9] where 
we observe that, in the initial stage, l c increases through a sequence of steps, each signaling the 
crossing of a corresponding barrier. Since the energy barrier decreases at increasing l c , the overall 
process speeds up as time unfolds. Summarizing, for £ large enough to guarantee that the equilib¬ 
rium probability distribution is given by (l4lb . the non linear effect due to (l20l) and the metastability 
of the region (j) —0, conspire to drive the system to a stable band compatible with a value of the 
stress o = <7 ( °)(£). Thus, upon increasing the shear, we should observe a linear decrease of the 
free-energy with the shear S (see equation (l28b ~). up to the point when the shear-band size l c be¬ 
comes of the order of the system size and the stable homogeneous solution is energetically favored 
again. This is again confirmed by the numerical simulations: in the top panel of Figure [10] we show 
the value of the free-energy for a Couette flow as a function of the apparent imposed shear S. We 
compare the space-averaged free-energy evaluated from the numerical simulations (red bullets), 
the free-energy for compact solutions (l28b (red dashed line), and the free-energy of the linear HB 
velocity profile ([2Tb (black solid line). The value of the noise amplitude is £ = 5 x 10 7 . In the 
inset of the same figure, we show the stress/shear relation obtained through the numerical simula¬ 
tions. A clear plateau at a ~ 1.5 is observed, corresponding to the value (J <0 ^ given in (l39l) . The 
velocity profiles are reported in the bottom panel of Figure [TOl 

Figure [TO] also demonstrates that the model defined by equations (l20l)- (l3TT) shows stable shear- 
bands for a < oy U \ where oy 1 is the value of stress in the upper plateau of the shear-stress rela¬ 
tion, i.e., basically cr' 01 in ([39b . In other words, the stability of shear-bands as a local minimum 
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of the functional F[(j)] and the noise amplitude £ cooperate in such a way as to increase the value 
of the yield stress at which (j) o becomes unstable. This cooperative effect highlights the subtle and 
somewhat counterintuitive role of the mechanical noise. 

Before closing this Section, a few comments on the role of q are in order. All of the above results 
have been obtained using the same value tq = 0.04 and it is natural to ask how the picture dis¬ 
cussed so far would change upon changing . According to (l39l) . one would naively expect that by 
increasing the value of o i{)] should also increase. However, a subtler analysis reveals that this is 
not the case. Since R ~ tq 2 ((d y ty) 2 ), there are two competing mechanisms which concur to fix the 
value of m (0 ' 1 ~ R 1 / 3 : the prefactor q 2 and the spatial average inherent to the definition of R. The 
former obviously increases with tq, the latter however has just the opposite effect. Indeed, since 
the system is split in two regions, the state (f>o and the compacton </>/.;, the spatial average involves a 
factor (1 — l c ) which is apparently decreasing upon increasing l c via an increase of q. As a result, 
the system develops a much weaker dependence on q than one would expect. 


VIII. DISCUSSION 


The formation of s 


ferent models 1119 . 


M 


lear-bands 


M 




las been investigated by several authors, using a variety of dif- 


5214651]. The permanent banding observed for viscosity-bifurcating 


YSF’s echoes the shear-banding effects that have been extensively studied in polymers and worm¬ 


like micelles [52 


58]. There, the criterion for the formation of steady state shear-bands hinges on a 


non-monotonic shear-stress constitutive relation, whereby the negative slope would trigger the in¬ 
stability leading to heterogeneous flows and shear-bands out of an underlying homogeneous shear. 
Mesoscopic models based on a non-monotonic constitutive relation have also been put forward to 


explain shear-bands in viscosity-bifurcating YSF’s 1119 . 


24, 


26 


61 


63], with the non-monotonicity 


that originates from different mechanisms. The well-known SGR (Soft Glassy Rheology) model 


proposed by Sollich and coworkers [ 161419] describes the jamming transition in terms of trap dy¬ 


namics for mesoscopic stress elements. Yield processes are taken as activated by some effective 


temperature x, based on the idea that yielc 


events elsewhere in the material cause “kicks” which 


add up to an effective thermal noise [16-180. The original SGR model captures a non-linear consti¬ 


tutive curve, yet a monotonic one, and thus fails to tie up with the shear-band instability mechanism 
discussed above. Shear-bands can however be recovered by lifting the assumption of a constant ef¬ 
fective temperature x, and considering a relaxation-diffusion dynamics for such parameter. For full 
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details see the original paper by Fielding et al. I119ll . An alternative model, by Mansard et al. [26Q, 


resorts to the KEP description [30Q, supplemented with a phenomenological equation, to account 
for the coupling between the flow and the structure. The authors show that viscosity bifurcation 


occurs, due to a flow-induced weakening of the structure. Simi 


and co-workers [66]. Another approach, due to Martens et. al. [611], uses a minimalist mesoscopic 


ar ideas, were studied by Dahmen 


model incorporating the local dynamics of plastic events (small-scale rearrangements, followed 
by stress redistribution). The authors studied the effects of the typical restructuring time needed to 
regain the original structure after a local rearrangement [630. The spontaneous formation of per¬ 
manent shear-bands is observed when this restructuring time is large compared to the typical stress 
release time. The formation of shear-bands due to structural relaxation has also been addressed 
in a detai 


ed mesoscopic study by Jagla [641 [65]. In the fluidity models considered by Picard 


et. al. [24], the authors study the relaxation dynamics of the fluidity via a free-energy functional 


approach, which includes ad-hoc deterministic terms to account for the competition between re¬ 
laxation (aging) and flow-induced rejuvenation. Such competition gives rise to stationary solutions 
as local minima (maxima) of the free-energy functional, thus leading to stable (unstable) branches. 
To be noted that the authors impose the stress as a spatially constant variable, thus suppressing the 
“mechanical” instability induced by inertial terms. Destabilization of states is however triggered 


by local stability analysis of the stationary free-energy solutions. Recent studies a 
shear-banding can arise quite generically under time-dependent flow protocols 167 , 


so show that 


68], even for 


materials with monotonic flow-curves. Once formed, these bands may persist only transiently, or 
may remain to a steady state, depending on whether the material’s flow-curve is monotonic or not. 


Further studies [140 . 


41 


69] point to the existence of a static yield stress exceeding the dynamical 


yield stress, with the latter measured upon decreasing the shear down to zero: this produces a step 
discontinuity in the material’s constitutive properties at zero shear and allows again the coexistence 
between unsheared and flowing bands. Yet another possible mechanism for genuinely steady state 
shear-banding without a non-monotonic constitutive curve, is the coupling between flow and con¬ 
centration [70]. Such mechanism has been explored in the context of steady shear-banded states 
of hard sphere colloidal glasses in Besseling et al. [71]]. 


The present work inscribes into the mesoscopic framework [16-28] and hinges on the idea of co 


operativity [29 h 32[] . All the results presented here stem from a free-energy picture (in analogy to 
equilibrium dynamics), a scenario that fits well within the framework of the simple models dis¬ 
cussed by Picard et. al. [24]. More precisely, the free-energy functional ([7]) has been designed 
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to lead (upon minimization) to the HB and fluidity results given in 11291. 


30], i.e., the “mean-field' 


scenario for the model. The built-in HB monotonic flow-curve does not admit permanent banding 
under conditions of a steady applied shear flow. However, the main point of the present work 
is that the addition of noise to the mean field scenario (|3TT) proves instrumental in promoting a 
regime of instability with negative slopes in the flow-curve for homogeneous flows (see Section 
m . We emphasize that also in other models described before [19, 61], permanent bands could 
be obtained only by extending the investigation beyond the mean field level. The introduction of 
fluctuations on top of the cooperativity scenario has the effect of decorrelating fluidity fluctua¬ 
tions on scales shorter than the cooperativity length %, provided the order parameter is sufficiently 
small. Thus, fluctuations decrease the average fluidity, a mechanism which competes with the non¬ 
locality triggered by cooperativity effects (non-local squared gradient terms in ([7])) and conspires 
to form stable bands. Our bands may be viewed as a “collection” of plastic events within a solid 


region that remains elastic (vacuum), thus echoing the works by Martens et al. 


6lh and Mansard 


et al. [260, without assuming neither a direct model of the restructuring time [61], nor an equation 
describing the flow-structure coupling in the fluid Qj. 

We also wish to emphasize the role played by diffusive terms. First of all, we remark that many 
rheological models have introduced diffusivity effects in the constitutive equations: in some mod¬ 
els (see [72] for a review) the diffusive terms are added to the stress equation and are needed to 
confer a finite interfacial width between the bands; other rheological models for shear-bands [60] 
propose an expansion of the stress tensor with the inclusion of second order derivatives of the 
shear, referred to as the “shear-curvature viscosity”. At variance with these models, we do not 
introduce any diffusivity effect in the stress, which is considered spatially constant [24]. Never¬ 
theless, we underline the crucial role played by diffusive effects in our model, which physically 
represent the spatial range of plastic relaxation. This effect is an essential mechanism to obtain 
heterogeneous compact solutions, as outlined in Section dill]) . Also, in our case, the diffusivity 

'feet of diffusion decreases and 


term is rather peculiar, see (1301) . In particular, for small / the e 


eventually vanishes at / = 0. This is not the case of other models [24]. We also emphasize that the 
negative slopes in the flow-curve result from the combined effect of both noise and cooperativity: 
failing either of these two, the renormalization of the flow-curve would be lost and no permanent 
shear-bandings would be observed. 


We finally wish to stress that our model falls in the general arena of mesoscopic models [27 


61 


730 


that do not include dilational effects which might take place during plastic events. These effects 
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may well add corrections to the picture drawn here, with the associated flow-concentration cou¬ 
pling phenomena discussed by Besseling [71fl. Concentration effects may possibly be introduced 
in the cooperativity length or in the yield stress [29[] of the model, but this would obviously imply 
the introduction of extra equations accounting for the concentration dynamics. This surely makes 
a very interesting topic for future work. 


IX. CONCLUSIONS AND OUTLOOK 


It is know that t 


be explained 


task 126, 


41 


7, 


le phenomenology of shear-bandings shows many intriguing effects to 
Modeling such a phenomenology is definitively a challenging 


11 

-14 

74] 

6L 

61- 

63, 


63-65 


67H69L 


7111 . In this paper we have explored the connection 


between shear-bandings and the cooperativity between local plastic events, characterized by a 
cooperativity length [ 29(43211. This fits in the framework of the fluidity models |29 h 32||. which 
have been proposed in the literature to explain the cooperative flow of complex fluids in confined 
systems. Technically, based on the simple observation that the fluidity is a non-negative definite 
order parameter, we have reformulated Bocquet el al. free-energy functional [30|] in terms of 
its square root (j) = if 1 / 2 , which is a signed quantity. For the geometry of a Couette flow, 
it is shown that once the stress-strain constraint Jq (j) 2 (y)dy — S/a is taken into account, this 
functional is minimized by inhomogeneous compact solutions (compactons), which coexist with 
regions of zero-fluidity (vacuum). Since the latter are unstable, the compactons increase their 
size at the expense of the vacuum, until they reach the size of the full system, at which point no 
further decrease of the energy can be achieved, other than by recovering the homogeneous HB 
solution, namely a linear Couette profile. This leads to a host of remarkable dynamical effects, 
primarily aging. Given the highly non-trivial dynamics of the fluidity field, we have further 
investigated how such dynamics is affected by the presence of stochastic fluctuations, typically in 
the form of mechanical noise. We wish to emphasize that such noise is inherently non-thermal 
in character and must be regarded as spontaneous fluidity fluctuations which occur also in the 
absence of any external load. To be noted that, since the fluidity is non-negative, so must be its 
average, which implies non-trivial restrictions on the noise strength. Such restrictions do not 
apply to the square-root fluidity, which, being signed, is free to fluctuate around zero. It is shown 
that, owing to non-trivial renormalization effects, if the noise is sufficiently strong, the unstable 
vacuum is stabilized, thereby paving the route to stable shear-band configurations. This qualitative 
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picture is confirmed by numerical simulations of steepest-descent dynamics under different initial 
conditions and noise amplitudes. The two starting ingredients of the present picture, i.e., fluidity 
as an order parameter and mechanical noise as a promoter of escapes from free-energy minima, 
are not new. The emerging picture, however, definitely is. In particular, it provides a transparent 
link between shear-bands and compactons, as well as a subtle stabilization mechanism via 
non-trivial noise renormalization effects. Finally, we remark that the above rich picture depends 
only on two free parameters: the cooperative length t, and the noise strength £. Many further 
directions for future research can be envisaged. As we have already pointed out earlier on, we 
have chosen a very simple form of the noise to compute analytically the renormalization effects 
of the free energy, and we believe the qualitative picture emerging is quite robust and independent 
(qualitatively) of the structure of the noise. This said, it is of course of decided interest to explore 
also the dynamic heterogeneity effects brought about by space/time correlations: recent work 
has highlighted strong correlations in the fluidity fluctuations [75r[77] and these correlations are 


likely to exert a significant effect on the flow response near the yielding point [77-81]. It is also 


of interest to consider the effects of shear-dependent fluidity fluctuations: note in fact that noise 
can also be triggered by the shear dynamics itself and through the elastic response of the material 
to local fluidity rearrangements. For example, one could enrich the noise with a dependency on 
the plastic activity, as proposed in the Hebraud-Lequeux framework [82] that underlies the KEP 
model presented in Bocquet et al. [30fl. 

Another important issue to be investigated in the present model is the role of boundary conditions. 
We wish to point out that the boundary conditions may serve as a trigger for the heterogene¬ 
ity [241], hence they have a strong bearing on the onset and stability of compact configurations. 
In the present work, we have focused on possibly the simplest boundary conditions compatible 
with compact solutions, although we expect that more general classes of compact fluidity 
configurations should arise in the presence of different boundary conditions. This is an interesting 
topic which warrants a separate study on its own for the future. 
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Appendix A: Stability of Compactons 


To simplify matters in the stability analysis we rescale the space and order parameter variables 


as 


so that the free-energy becomes 



0 = m(j) 


(Al) 


F[$] 





dy. 


Then equations (fl3l) and ([131) become, respectively: 


(A2) 


j 2(j) 2 dyy(j) + 2 <j) (<5y0) 2 + 0 3 - <j ) 3 0 | = 0 

|<p(%0) 2 =£-±0 4 + i<F 


where the constant E has been rescaled by a factor m , i.e., E = E/m 5 . To assess the stability of 
the compactons fa(y), we compute 8F = F[fa + <50] — F[fa\ up to second order terms in 8(j). 
Based on (l24l) . stability is guaranteed if 8F is positive defined. After a rather straightforward 
computation we obtain (fa > 0) 


0/ ~ ryo+2lc _ „ 

8F = 2m 9 / 2 ^ / (8<t>Y 

Jy o 


(dyfc ) 2 ~ ^$e + $E 


dy + F- 


+ 


(A4) 


where F+ includes all positive defined terms 


F + = 2m 9/2 q f 

Jy, 


yo+2l c 


(j> 2 (dy8(j)) 2 dy. 


y o 


The integral in (IA4I) is computed in the interval [yo,yo + 2 l c \ where the localized solution is de¬ 
fined. Because of symmetry, it is enough to show that the integral between [vq.vo + l c ] is positive 
defined. Furthermore, we can always choose the origin of integration with the position yo where 
the localized solution fa becomes zero and write 


kiy) = m'/«A(y)y '/ 2 


where A is an analytic function of y. Since AdyA -\-y(dyA ) 2 > 0 for 0 <y<i c and F + is positive 
defined, we obtain 


8F 

4 m 9 / 2 <§ 



■+f' 2 (AE) l / A A(y) 


dy. 
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We have then studied the properties of A(y), finding a lower bound for it: A(y) > 1/ y/2. Hence 




±_Z + Ai(4£)l/4 
4y 2 + V2 l ' 


dy. 


(A5) 


Because the function [A — 1 +y 3 / 2 (4£ , ) i / 4 -fj] is positive defined for E >E C = 0.01, it follows that 
all localized solutions with integration constant E in the range E e [0.01,0.05] are local minima of 
F[(j>], where E = 0.05 corresponds to the upper bound Em = m 5 /20 discussed in SectionUni above 
which no compact solution can be found. It is clear from our result that the number of possible 
stable compactons is still very large. 

Appendix B: Self-Consistent Field approximation 


In this appendix, we derive the expressions (f39l)-(l40b discussed in Section (Uni) . We start with 
equation OTI) 

0 = — 

F[^} = 2 fo [-?m0 4 + i0 4 |0| +% 2 (/) 2 (dy < l>) 2 ]dy 

which we rewrite in the form 


(Bl) 


M = -2^+4^ 2 0 2 ^+4^ 2 0(^0) 2 + v / ew (y,t) 


(B2) 


4 r ' 5 ' 

Our aim is to show that, because of the noise, V is renormalized and, moreover, the renormalization 
is state dependent. The difficult terms to be estimated are those proportional to t, 2 in equation 
(IB2I) . i.e., (j) 2 r) vv 0 and (j)(d y (j> ) 2 . Starting from the free-energy, a self-consistent Hartree-like field 
approximation [49451] can be performed on the term (t> 2 (d y <f >) 2 


F[0] = 2 


-h"^ 4 + j^ 4 l*l+f * 2 + j(M ) 2 


dy 


(B3) 


where 


R = 2^({d y ^) 1 ) D = 2^((j) z ). (B4) 


Note that formally we should apply the same approximation to the terms (j ) 4 and 0 4 \d\ but for our 
purpose this would add only small perturbations to the results below described. Next, by using 
equation (IB3I) . we can rewrite equation (IB2I) as 


— —2—^- +2 Ddyy(j) + y / £w(y,t) 


(B5) 
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where the effective potential V eff is given by 


Kff(0) = -jM</) 4 + ^</> 4 l</>\ + ^</> 2 . 


(B6) 


According to (IB6I) . the effective potential has local minima which can be found by solving the 
equation (0 > 0): 


dV, ff 

dif) 

d 2 V efi 

df- 


— (R0 — m0 3 + 0 4 ) 

4> = 4> m ,n V ’ 

> 0 . 

0 ~0min 



0 


(B7) 


It is easy to show that (IB7I) has one solution 0 mm — 0o — 0 for m < m c ~ (21R/A ) l,/3 , whereas for 
m > m c there exist two local minima at 0 mm — 0 q — 0 and (j) mm — 0 m ~ m — R/m 2 + 0((R/nr) 1 ) (see 
top panel in Figure [5]for a sketch). To close the problem, we need to compute R = 2q 2 ((d y 0) 2 ). 
This computation can be done perturbatively (i.e., assuming £ to be small) by linearizing equation 
(1B6I) near the minima. 

We start with the minima at 0 — 0,„ ~ m — R/m 2 + 0((R/m 2 ) 2 ). As we shall see in following, R 
is proportional to £ and we can simplify the computation by assuming 0 — m. Next, we consider 
the fluctuations 80 near 0 — m which obey the stochastic differential equation 


d t 80 = 4 t > 2 m 2 d yy 80 — 2(m 3 +R)80 + ffe w(y,t). 


(B8) 


The above linear equation can be solved in Fourier Transform for 80k and we obtain 

d t 80k — 4^ 2 m 2 k 2 — 2m 3 — 2R)80k+ \/ew(k,t) 


(B9) 


where now both 80k = —= / exp (iky)80(y,t)dy and w (k,t) are complex variables and moreover 
(w(ki,?i)w*(k 2 T 2 )) = 8{k\ +k,2)8{t\ —ti). Equation (1B9I) is a linear Langevin equation for each 
k whose asymptotic variance can be easily computed. After a simple algebra we obtain: 


R = 2 % 2 f k 2 {80 k 80l)dk=-^ 

ek M 


t; 2 k 2 


m 


2b, 2 k 2 + m + R/m 2 
2 (1 + 0(m/kM) + 0(k M R/m 2 )). 


dk 


Upon recalling that m/. ~ R and substituting for R the result obtained in (IB 101) . we obtain m/ 


(BIO) 


ekM 


and thus 

m c ~ (£k M ) 1/5 . (Bll) 

We refer to this value of m c as nff ] . Our derivation implies that for m < m2 1 ' the only minimum 
of the effective potential is 0 q while for m > there exist two local minima, and in particular 
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the minimum near 0 = m is the global minimum of the potential. 

Next we consider the minima 0 = 0. In this case the problem is more complicated since R should 
be computed self-consistently with D. To solve the problem, we linearize equation (1B6I) near (j> — 0 
obtaining: 

^50=2Do^50-2/?50 + Vew(y,f). (B12) 


We next perform a Fourier transform obtaining a set of linear Langevin equations. A simple 
algebra gives: 

D = 2$ 1 J(S<l> i 5r k )dk=2sl; 2 J (B13) 

r f k 2 dk 

R = 2$ 2 J ^(S^S^dk = 2e % 2 I Dk 2 +R - (B14 ) 

By multiplying equation (IB 13b by R and equation (IB 141) by D, and summing the two results, we 


get 


2DR = 2eE, 2 k M . 


(B15) 


Another relation between D and R can be explicitly obtained by solving exactly the integral in 
(IB 13b and assuming k M large 


D = 2e^ 2 ! 
£<§ 2 


dk 


= 2 - 



= 2 


~R 


dk 


R 


-atan 


Dk 2 /R + 1 
7T£<§ 2 [R 



(B16) 


The self-consistent solutions of (IB 1 5b - (IB 161) are given by: 

D ~ ^ M ; R~ %kM \J£kM■ 

k-M 


(B17) 


Equation (IB 171) implies that the fluctuations near (f )o are given by 

<W)2> = ?~W (B18) 

meaning that near <j )o the average fluidity / 0 is not zero, which is the result reported in (1331) . Finally, 
we recall that m 2 ~ R, so that 

m c ~ (^mv / Sm)' /3 (B19) 


and we refer to this value of m c as m {{]] . 
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In order to gain a physical insight in the above approximations, it is expedient to consider again 
equation (1B2I) and decompose (f) as follows: 


0 = ( 0 ) + ^ 


(B20) 


where 8(j) are (small) fluctuations around the space average ((f)). Equation (IB2I) is then rewritten 
as: 

= 2/n^ 3 -20 3 |0|+4^ 2 ^ > ,[0 2 (<5 v ^)] -4<§ 2 0(<3v^) 2 +v / ew(y,t) (B21) 

Upon performing the space average of (1B21I) and assuming self-averaging (/.<?., ensemble average 
is equal to space average), we obtain: 

(B22) 

where we have neglected terms proportional to (8<f>)(d y 8(t)) 2 and also terms coming from the 
decomposition of (j) 3 and 0 3 |0|, which are supposedly much smaller than (d y 8(j)) 2 . By linearizing 
the rhs of (1B21I) around ((j)) >0, and subtracting equation (IB221) . we obtain (to the linear order in 
8(j)) the following equation 


d t 8<$> = 2(3m(0) 2 -4(0) 3 )<5<)> +4<§ 2 (0) 2 <9yy<50 + v / ew(y,t). 


(B23) 


Next, assuming that (0) ~ m + 0(e), we can compute the quantity R — 2^ 2 ((d y 8(j)) 2 ) from (IB231) : 
the equation is linear and can be solved similarly to what we have done for (IB8I) . obtaining the 
result in (IB 101) . 

Moreover, we can recover the results for by writing equation (1B2I) for the fluidity f — (j) 2 . 
Upon using Ito calculus [83], we obtain: 


d t f — £k M + 0(f 2 ) + t, 2 fd yy f + y/i/w(y,t) (B24) 

where the terms 0(f 2 ) are neglected in the following since we are investigating the stability of the 
state (j) — 0. Upon averaging in space (IB24I) and assuming self-averaging to hold, we find 


dt(f) = -Z 2 ((dyf) 2 ) + £k M . (B25) 

We then use the approximation ((d y f) 2 ) — ak 2 M (f) 2 , where a is a constant of order 1. At the 
stationary state, equation (1B25I) predicts 


$ 2 ak 2 M (f) 2 = ek M 


(B26) 
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which is equivalent to 


<W) 2 > = (B27) 

CC l '~C,kM 

Equation (IB27I) . apart for a numerical constant, gives exactly the same result reported in (IB 171) - 
(lBl8l) . 
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FIG. 1: Properties of compact solutions <j) (: (fl4l) . Top panel: a sketch of a compact solution obtained from 
equation (fl3l) : the order parameter 0 is non zero only in some compact sub-domain [yo. yi where the 
profile is identified with the compacton <f>£ (red dashed). The compacton </y.- has a maximum <j)(y) = <j) 
in v = (yo +y\)/2. Outside the compact sub-domain the order parameter is identically zero (blue dashed- 
dotted). Bottom panel: the average fluidity associated with the compacton i.e., f ( ] <p/(y) dy, for different 
values of the imposed stress a in the free-energy (fl2l) . The average fluidity is reported as a function of the 
constant of integration E (see equation (fl5l) ) normalized to its maximum allowed value Em(<j) = ^ 

(see equation ( fl6l) and text for details). Notice that the average fluidity is made dimensionless with respect 
to the inhomogeneity parameter q of the free-energy (fl2l) . In the inset we report the half width l c of the 
compacton calculated according to equation (flSl) as a function of E/Em, again for different values of the 
imposed stress a (same as the main panel). 
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FIG. 2: Top panel: we show the global free-energy F[(j)] (see equation (1121 )) for three different cases: linear 
velocity profile corresponding to a HB law (black solid line); compacton adjacent to the upper wall of 
the channel (see also bottom panel) with stress a = 1.4 (blue dot-dashed line) and a = 1.5 (red dashed 
line). At a given stress a, the size of the compacton is changed in order to satisfy the constraint (l20l) . The 
free-energy of the compactons is smaller than the linear velocity profile up to a critical shear where the 
size of the compacton becomes equal to the channel width. Bottom panel: free-energy densities L((j),dy(j)) 
defined in equation (fl2l) as a function of the channel position at imposed shear S = 0.1 (subscripts indicate 
the compact (c) and the Herschel-Bulkley (HB) solutions). The chosen compact solution (red solid line) 
corresponds to a stress a = 1.5. In the inset we report the corresponding velocity profiles, for both compact 
and HB solutions. 
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FIG. 3: The behavior of a(t,t w ) under the effect of the shear S = 0.1 applied to the system after the waiting 
time t w . The initial condition for 0 are: 0 = 0.2 for y € [0 : 0.8] and (j) = 0.5 for y € [0.8,1], The black 
arrow indicates the time 0,(t w ) at which the value of the stress reaches the HB value 1 + S 1 / 2 . In the inset we 
show ti as a function of t w : a clear logarithmic dependence is observed. Details of numerical simulations 
are reported in the text (see Section [Vj. 
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FIG. 4: We report the average fluidity /o = ((j) 2 ') (see also equation (l33l) i obtained from numerical simula¬ 
tions of the model equations (I3TT) . Main panel: the average fluidity as a function of q for £ = 10 8 . Inset: 
/o as a function of the noise strength £ for q = 0.04. In both cases the scaling predictions of equation (1331) 
are verified. All simulations have been performed at constant stress a = 1 (i.e., m = 0 in (l3Tlfl . Details of 
numerical simulations are reported in the text (see Section Ivl). 
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FIG. 5: Top panel: The effective potential V cll defined in equations (l36l) with R = 0.001. A bifurcation 
is present at m = m* = (27R/4) 1 ^, above which two local minima appear at 0 min = 0o = 0 and at (j) mm = 
0,„ ~ m — R/m 2 + 0((R/m 2 ) 2 ), together with a local maximum in 0 = 0m- Note that there is a range of m 
(m* < m < m c ) where the local minimum 0 m appears, but 0o is still the global minimum. At the critical value 
m = m c , the vacuum and homogeneous (stable) fluidized solution can coexist. Bottom Panel: homogeneous 
stationary flow-curve obtained from the local extrema equation R — ra0 2 + 0 3 = 0 (see text for details). Due 
to fluidity fluctuations (R > 0), such flow-curve is a decreasing function up to the point in which it slopes 
upward at higher shear. This indicates the onset of an instability (see text for details). 
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FIG. 6: Top panel: Stress shear relation obtained by the numerical simulations of (I3TI) with £ = 10 H . 
Starting with an homogeneous initial condition (j) = (j) 0 =Owe increase - step by step - the stress (red 
circles) up to the value where a clear HB behavior (solid line) is detected. Next we decrease the stress 
starting with 0 = 0,„ = m (blue circles). A clear hysteresis cycle is observed. Arrows indicate the points 
where the values of and m^^(e) are extracted. Correspondingly, the critical values of the stresses 

can be obtained based on (I39l)- ll40l) . Bottom panel: same as the top panel with £ = 5 x 10 1 . 
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FIG. 7: We report the values of mS Q \e) (red circles) and m 111 (e) (blue squares) computed from the numer¬ 
ical integrations of (|3TT) with noise amplitude in the range [ 10 8 : 5 x 10 7 ]. The numerical procedure to 
determine the transition values, m^(e) and m (1) (e) , is highlighted in Figure [6] 
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FIG. 8: Top panel: time dynamics of a pre-sheared state under the influence of a constant stress a based on 
equation (f3Tb . We consider an initial homogeneous state </>(y. 0) = (j) a corresponding to a pre-sheared state 
with shear 5(0) = S a = ty 2 o = 0.01. We then study the behavior of the apparent shear 5(f) = <7 /,, (p 2 (y,t)dy 
for different values of the applied stress a. There is a clear bifurcation between the values of the stress 
below a critical value a^ c \ at which the shear tends to vanish, and those above o ! ' c \ where the shear goes 
asymptotically to a non zero value. Bottom panel: we show 5(f) for different values of the initial shear 5„ at 
fixed applied stress a = 1.2. For small values of S a , the apparent shear tends to vanish whereas for S a greater 
than some critical value, the shear goes asymptotically to a non zero value. The numerical simulations are 
based on the model equations ( |3TI) with e = 10~ 8 , £ = 0.04. 
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FIG. 9: Top panel: we show the band size l c (t) as a function of time in a Couette flow with two different 
initial conditions, l c ( 0) = 0.20 and / c (0) = 0.75, and 2 different noise strengths, e = 10 7 and e = 10 s . 
Bottom panel: we report the band size l c (t) as a function of time with initial conditions chosen in the range 
[0.18: 0.22] and noise strength e = 10 H . The inset reports the early stages of the time dynamics for the case 
with l c ( 0) = 0.18, showing the increase of the band through a sequence of steps. In all cases the imposed 
shear is S = 0.04. 
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FIG. 10: Top panel: we show the free-energy for a Couette flow as a function of the apparent imposed shear 
S and for £ = 5 x 10 1 . The red bullets correspond to the space-averaged free-energy F[(j)] in the stationary 
state obtained from numerical simulations of the model equations (OTI) : the red dashed line is the value of 
F for compact solutions (1281) : the black solid line is the value of F corresponding to the linear HB velocity 
profile (OTI) . In the inset we show the stress/shear relation obtained by the numerical simulations. A clear 
plateau at a « 1.5 is observed corresponding to the value cr l0; of equation (l39l) . Bottom panel: the velocity 
profiles corresponding to the shear-band region of the top panel are reported. 
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